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THE SPIN HISTORY OF PROTOSTARS: DISK LOCKING, REVISITED 
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ABSTRACT 

In this talk, we take a new look at the theory of disk locking, which assumes that an accreting protostar rids 
itself of accreted angular momentum through a magnetic coupling with the accretion disk. We consider that 
differential rotation between the star and disk twists the field lines. For large enough twist, the magnetic field 
lines connecting the star and disk open and disconnect. This significantly reduces the spin-down torque on 
the star by the disk, and so we find that disk-locking theory predicts spin periods that are much too short to 
account for typical observed systems. 
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1. INTRODUCTION 

The collapse of a molecular cloud naturally leads 
to a phase consisting of a central protostar sur- 
rounded by a centrifug ally supported accretion disk 
(for a review, see lBodenheimerlll995|) . Disk winds, 
MRI turbulence, and/or viscous processes remove 
angular momentum from the disk and results in the 
accretion of material with high specific angular mo- 
mentum onto the star. For typical parameters for ac- 
creting protostars (CTTS's), the accretion alone will 
spin the star up to near breakup speed in less than ~ 
10 5 years, assuming the star hadn't already formed 
at near breakup speed. Since the accretion lifetime is 
often greater than 10 6 yr, the stars must rid them- 
selves of this excess angular momentum. Further, 
it has been generally accepted that accreting proto- 
stars spin mor e slowly than their n on-accreting coun- 
terparts (e.g.. iBouvier et aTlll993lh The general ex- 
planation is that the presence of an accretion disk 
somehow regulates the stellar spin, and then after 
the disk is dispersed, the star spins up as it contracts 
to ward the main sequence. 

Konigl (199l |) applied the neutro n star accre- 
tion model of iGhosh fc Lambl (|l979|) to accreting 
protostars and showed that a "disk-locking" (DL) 
mechanism could explain the coincidence of accre- 
tion and slow rotation, in those systems. Ac- 
cording to DL theory, magnetic field lines connect 
the star to the disk (acting as "lever arms") and 
carry torques that oppose and balance the angu- 
lar momentum deposited by accretion. CTTS's are 
now known to posses kilogauss-strength fields (e.g., 
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Fig. 1. Magnetic star-disk interaction. (from 

iMatt fc Pudritj|2004h 



x McMaster University, Canada. 
2 CITA National Fellow. 



Ijohns-Krull. Valenti. fc Koreskoll999|) . and the gen- 
eral DL model has been invoked by many authors. 

Recent ob s ervati ons of CTTS's in Orion by 
IStassun et alJ l)l999|) . however, show no correla- 
tion between observed rotation periods and accre- 
tion diagnostics, calling the standard DL scenario 
into question. Furthermore, the magn etic fields of 
CTTS's, while strong are disordered l|Safierlll998l 
I.Tohns-Krull et al . 1999), which reduces the effective- 
ness of magnetic torques required for DL. These de- 
velopments prompted us to revisit the general the- 
ory of DL. In particular, the connectivity between 
the star and disk is an important issue. Much re- 
cent work has shown that the magnetic connection 
between the star and disk is severed when the mag- 
netic field is highly twisted. Here, we show that the 
resulting spin-down torque is significantly reduced, 
and the DL model cannot account for accreting stars 
that spin slowly (e.g., ~10% of breakup speed). 
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2. THE "STANDARD" DL MODEL 

To begin, we first formulate a basic model that 
is representative of the general DL picture dis- 
cussed by many authors. We follow the work of 
lArmitaee fe Clarke! l|l996t hereafter AC96). The 
general theory assumes the central star contains an 
axis-aligned dipolar magnetic field. A dipole is re- 
quired because the field strength falls off as slowly 
as nature allows {r~ 3 ), and any higher order mul- 
tipole falls off so quickly that torques become neg- 
ligible. The dipole field is anchored in the surface 
of the star and also connects to the accretion disk, 
which is assumed to always be in Keplerian rotation. 
The disk accretion rate M a is constant in time and 
at all radii in the disk. Figure illustrates the ba- 
sic idea and identifies the location where the disk 
is truncated (Rt), the outermost radius where the 
closed stellar field is connected to the disk (iZ ut), 
and the corotation radius (R co ), where the Keple- 
rian angular rotation rate equals that of the star. 
The usual assumption is that R out ^> R co (AC96 
used Rout —> oc). 

The magnetic torque on the star from field lines 
threading some range of radii in the disk midplane 
is given by 

r ln = / l—r^r where 7 = — - (1) 
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(e.g., AC96) where \i is the dipole moment, and 7 
is the "twist" of the magnetic field. So the torque 
depends not only on the strength of the field but, 
more importantly, on how much it is twisted. The 
more it is twisted (larger 7), the stronger the torques. 

The field twist is generated by the differential ro- 
tation between the star and disk. As the field is 
twisted, it resists the twisting (hence the torques) 
and slips backwards through the disk. The larger 
the 7, the faster the slipping. When it can slip back- 
ward at the same rate as the differential rotation 
rate between the star and disk, a steady-state for 7 
is achieved. The speed of slipping field lines is given 
by Vd — /3t>kcp7, where (3 is a constant scale factor 
by which compares to the Keplerian speed, Wkop 
(AC96 use = 1). Thus, the steady-state configu- 
ration of 7(r) is given by 7 = 0~ 1 [(r/R co ) 3 / 2 — 1], 
and so the torque in equation ^ can be calculated. 
The value of is unknown . Standard a-disk physics 
jShakura fc Sunvaevl fl973^ gives an upper limit of 
< 1 and a likely value of a few orders of mag- 
nitude lower. We consider a value of = 10~ 2 as 
reasonable, but given the uncertainties, we keep 
as a free parameter in our analysis. At first, we will 
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Fig. 2. Differential torques in the disk midplane for a 
system with /3 = 1, M a = 5 X 10" 8 M Q yr _1 , M» = M Q , 
R, = 3R Q , B, = 10 3 G, and a stellar spin period o f 5.7 
days (so R co = 4.5R*). ffrom lMatt fc PudritzlEuOl 

use 0=1 (which gives the solution of AC96), but 
we consider other values in the next section. 

Figure|5]shows the differential torques (per Sr) as 
a function of radius (normalized to R co ) in the disk 
midplane, for a system with the parameters listed in 
the figure caption. The line labelled "accretion" rep- 
resents torque that is required to supply the steady- 
state M a . The line labeled "magnetic" shows the 
differential torque (from eq. from the stellar field 
threading the disk. Inside R co , this torque acts to 
spin the star up. It decreases rapidly away from the 
star as the dipole field becomes weaker. It is zero 
at R co because the differential rotation (and thus 
7) is zero there. Outside R co it becomes stronger 
(though now spinning down the star) as 7 increases, 
but it eventually becomes weaker again because the 
decrease in the dipole field strength decreases faster 
than the increase of the 7. In order to satisfy the 
steady-state condition, the disk must restructure it- 
self so that the sum of the magnetic and internal disk 
differential torques must equal the accretion differ- 
ential torque at all radii. The dashed line labelled 
"viscous" in Figure|21shows the required internal disk 
torque. 

The disk will be truncated near where the accre- 
tion and magnetic differential torques are equal (and 
where viscous torque = 0). From that point (Rt) in- 
ward, all of the specific angular momentum of the 
disk material will end up on the star. So, to calcu- 
late the net torque on the star from the accretion 
of disk material, r a , one integrates the differential 
accretion torque (shown in Fig. |2J) from R t to the 
surface of the star. The net magnetic torque, r m , is 
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obtained by integrating equation ^ from R t to i? ut 
(which is thus far assumed to be oo). 

For any given values of M*, i?*, B*, M a , and the 
stellar rotation period, this "standard" theory gives 
the net torque on the star. The system is stable 
in that, for fast rotation, the net torque spins the 
star down, and for slow rotation, the star spins up. 
Also, for typical CTTS parameters, the torques spin 
the star up or down in ~ 10 5 yr, so one expects 
that most systems will exist in a spin equilibrium 
state where the net torque on the star is zero. The 
"standard" DL model thus predicts the spin period 
in the equilibrium state, T cq , at which the system is 
"disk locked." Figure [3 is shown in its equilibrium 
spin state (T eq = 5.7 d). 

Models such as this have had success at explain- 
ing the spin rates of slow rotators. For example, the 
well-studied CTTS, BP Tau has M = 3 x 1Q- 8 M M 
yr" 1 , R* = 2R Q , and M* = O.5M ijGidlbring et all 



Il998i). Using the mean field strength of 2 kG found 
bv lJohns-Krull et al.l (fl999). our "standard" DL the- 



ory predicts T eq = 7.5 d (corresponding to 6% of 
breakup spee d) — remarkably s imilar to the observed 
value of 7.6 d l|Vrba et alJl986h . Thus the DL theory 
seems to work, but there's at least one major prob- 
lem. Namely, we assumed that the star and disk 
were connected to i? ou t — > oo. At large radii, the 
field will be highly twisted, and there is a physical 
limit to that twist, which we have so far ignored. In 
the next section, we consider the effect on the DL 
model of an upper limit to the magnetic twist. 

3. EFFECT OF LIMITED TWIST 
Twisted magnetic fields that connect the star to 
the disk exert torques between the two. The larger 
the region in the disk that is magnetically connected 
to the star, the larger is the total magnetic torque. 
So the actual location of i? ut is important, since 
it delimits the connected region (i.e., it determines 
the int egration range in eq . [T|l . Many re cent stud- 
ies fsee lllzdenskv. Konigl. fc Litwinll2002l and refer- 
ences therein) have shown that dipole magnetic field 
lines connected to a rotating disk open up, when 
twisted past a critical value of 7 = j c « 1. This is 
unavoidable and occurs because the magnetic pres- 
sure associated with the azimuthal component of 
magnetic field pushes out against the poloidal field 
lines. These open field lines (see, e.g., the field lines 
outside R out in Fig. ^) cannot convey torques be- 
tween star and disk. 

Since the steady-state value of 7 increases as r 3 / 2 
away from R co (see it inevitably reaches the 
critical value, which we take as 7 C = 1. As an ap- 
proximation, we assume that the star is connected 




Fig. 3. Logarithm of the net magnetic torque as a func- 
tion of log(/3). The torque (denoted Xi see text) is nor- 
malized to the AC96 value, so that the value of x = 1 is 
the torque from the "standard" DL theory, and the solid 
line shows the effect of limited magnetic twist. For ref- 
erence, the dotted lin es show x = P an d X = P~ ■ (from 
iMatt fc Pudrit j|2004l) 

as before to the disk, but that the extent of con- 
nected region is now limited to the region where 
7 < 7 C . The outermost location of this region is 
-Rout = (1 + Plc) 2 ^ 3 Rco, beyond which, we assume 
the star is disconnected from the disk, and so the 
differential torque is zero. Since the size of the 
connected region is smaller than for the "standard" 
model, the net magnetic torque is less. 

Figure 03 illustrates how the magnetic torque is 
altered by the opening of field lines, as a function 
of the parameter (3. Shown is the value of the net 
magnetic torque (denoted x)i renormalized so that 
the torque predicted by the "standard" DL theory 
of the previous section {[3 — 1 and R out = 00) gives 
a value of x — 1- The normalization allows Figure |3 
to be valid for any given values of M a , M*, /i, and 
the stellar spin period. It is evident that the mag- 
netic torque will always (for any [3) be significantly 
less than predicted by the "standard" theory. The 
dependence of x on P can be understood as a com- 
petition between two different effects: One is that 
i? out decreases for decreasing (3, reducing the inte- 
gration range of equation ^ the other is that the 
steady-state 7 decreases for increasing (3, reducing 
the differential torque at each radius. For the criti- 
cal value of (3 = 1 , these two effects conspire to give 
a maximal value of r m that is four times less than 
predicted by the "standard" DL theory of section [21 
(for any given values of M a , M*, f2*, and fi). So by 
using (3 = 1 above, we have considered the "best pos- 
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sible case" for DL theory, since r m is less for all other 
values of (3. A reduced magnetic torque, means that 
the star must spin faster before it is in equilibrium. 
A faster spin reduces R co , so the torques come from 
closer to the star where the dipole field is stronger, 
making up for the decreased integration range. 

We can now revisit our example case of BP Tau. 
Using the "updated" DL theory (7 C = 1), the "best 
case" value of (3 = 1 predicts T cq = 4.1 d. For (3 — 
0.1 (or = 10), T cq = 2.5 d. The time to spin up 
from 7.6 d to 4.1 d (or even to 2.5 d) is 1 x 10 5 yr, 
which is signifi cantly shorter th an BP Tau's age of 
6 x 10 5 yr l|Gullbring et alll998lb Therefore, BP Tau 
has either gone through a recent change in (e.g.) M a , 
so that it is not currently in the equilibrium state, 
or the model is incomplete. In order for BP Tau 
to currently be in an equilibrium spin state, there 
must be significant spin-down torques on the star 
other than the torques along field lines connecting it 
to the disk. As an aside, we note that we have thus 
far used the mean stellar magnetic fie ld strength of 2 
kG found bv lJohns-Krull et all Jl999), though this is 
not the true strength of the dipole field. If instead we 
use the 3a upper limit to the d ipole field strength of 
200 G l| Johns-Krull et al.ll999h . even the "standard" 
theory predicts T eq — 1.0 d. 

4. SUMMARY OF PROBLEMS WITH DL 

The disk- locking scenario has recently been called 
into question by observations, as well as by theo- 
retical considerations. In particular, there are four, 
completely independent issues: 

1. IStassun et alJ {l999) found no correlation be- 
tween accretion parameters and spin rates of 
TTS in Orion. 

2. CTTS's a pparently do not have strong dipol e 
fields re. g ..ISafierll99stl.Iohns-Krull et a,lll999lb 

3. Stellar winds are expected to open field lines 
that would otherwise connect to the disk ijSafierl 
1998). A disk wind would have a similar effect. 

4. Finally, a large portion of the magnetic field con- 
necting the star to the disk will open up, due to 
the differential rotati on between the two (e.g., 
llJzdenskv et alJf2002|) . We have shown that the 
resulting spin-down torque on the star by the 
disk is less (by at least a factor of four and pos- 
sibly by orders of magnitude) than calculated 
by previous authors. The predicted equilibrium 
spin rate is therefore much faster. 



So the DL scenario does not explain the angular 
momentum loss of the so called slow rotators — the 
group originally targeted by DL theory. 

We conclude that, in order for accreting proto- 
stars to spin as slowly as 10% of breakup speed, there 
must be spin-down torques acting on the star other 
than those carried by magnetic field lines connecting 
the star to the disk. The presence of open stellar 
field lines naturally leads to the possibility that ex- 
cess angular momentum is carried by a stellar wind 
along those open lines. We plan to investigate this 
possibility in the near future. 

We'd like to thank the organizers for a fantastic 
meeting. This research was supported by NSERC, 
McMaster University, and CITA, through a CITA 
National Fellowship to S. M. 
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